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Abstract 

A novel method, connecting the space of solutions of a linear differential equa- 
tion, of arbitrary order, to the space of monomials, is used for exploring the 
algebraic structure of the solution space. Apart from yielding new expressions 
for the solutions of the known differential equations, the procedure enables 
one to derive various properties of the orthogonal polynomials and functions, 
in a unified manner. The method of generalization of the present approach 
to the multi-variate case is pointed out and also its connection with the well- 
known factorization technique. It is shown that, the generating functions and 
Rodriguez formulae emerge naturally in this method. 



*gurappa@ipno. in2p3.fr 
^panisp@uohyd.ernet.in 
•'■panisprsQuohyd.ernet.in 



1 



I. INTRODUCTION 



It has been recently shown, by two of the present authorsB, that the solution of a linear 



differential equation (DE), of an arbitrary order, can be mapped to the space of monomials, 
if the operators, relevant to a given DE, can be separated into a part containing the Euler 
operator {D = xd/dx) and the constants and another one retaining the other operators. 
Separating a given DE into two different parts and generating a series solution, from the 



advantage of the present method derives from the fact that, the Euler operator is diagonal in 
the space of monomials and any other differential operator or monomial (O'') is characterized 
by a definite degree d, with respect to the Euler operator i.e., [D, O'^] = dO'^. By a judicious 
use of these two results, the space of solutions of the above DE of arbitrary order, is directly 
connected to the space of monomials, avoiding any explicit integration, unlike the previous 
approaches. This technique, not only yields novel expressions for both polynomial and 
functional solutions of the known DEs, but can also be straightforwardly extended to a 
wide class of multi-variate DEs. It can be used to find solutions of equationsi, involving 
the so called Dunkl derivativesi; these type of equations are being extensively studied in 
the current physics and mathematics literaturei. Furthermore, the approach leads to the 
diagonalization of various many-body interacting systems0'ii of the Calogero-Sutherland 



In the present paper, we confine ourselves, to the single variable case and explore the 
origin of various algebraic structures in the space of solutions of the DEs and only briefly 
outline the procedure to generalize this approach to the many-variable cases. These alge- 
braic structures manifest transparently here, since they appear naturally in the space of 
monomials and our approach connects these monomials to the solution space. It is worth 
mentioning that these algebras are responsible for symmetries and various degeneracies of 
physical problems!, associated with the DE under consideration. 

The paper is organized as follows. In the subsequent section, we briefly outline the 



solution of the simpler one. 



by integration has been tried earlier in the literatureH. The 
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essential steps of our method and use the same to obtain novel expressions for the solutions 
of the well-known hypergeometric and confluent hypergeometric equations and also show 
its applicability to generalized hypergeometric equations. The ladder operators and the 
underlying algebras are then obtained in section III. We then establish the connection of 
this approach with the factorization technique and Supersymmetric quantum mechanics 
(SUSY-QM). The utility of the method, in finding various properties of the orthogonal 
polynomials, is then demonstrated in section IV, by finding the generating functions and 
the Rodriguez formulae, for a number of cases. We conclude in section V, after pointing out 
a number of problems, where the present method can be profitably employed. 

II. MAPPING BETWEEN THE SOLUTION SPACE OF DIFFERENTIAL 

EQUATIONS AND MONOMIALS 

In this section, we reproduce, for the sake of completeness, a recently proposed method 
of solving linear differential equations of arbitrary orderS and explicate the same, with a 
number of examples. 

After suitable manipulations, if a DE can be cast in the form 



[F{D) + P{x,d/dx)]y{x) = 



(1) 



where, D = x^, F{D) = Y,n=^oo ^nD'^ , o„'s are some parameters and P{x,d/dx) is a 



function of x, d/dx and other operators, then the solution to the DE can be written as, 

1 



yix) = cAY: (-1) 
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(2) 



provided, F{D)x^ = 0. Here Cx is a constant. The proof is straightforward and follows by 
direct substitution: 
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(3) 



Eq. connects the solution of a given DE to the space of monomials. It needs to be 
emphasized that, the inverse of F{D) is well defined in the above expression, since it is 
diagonal in the space of monomials. 

We illustrate the working of this method in the context of the well-known hypergeometric 
differential equation (HGDE)§, given by, 

d^ , ^ d ^ d^ d 



X + (a + /3 + 1) X— + a[3 - x-—t - 7— 



F{a,(3;r,x) = 



(4) 



dx"^ dx dx"^ ' dx 

Collecting the part containing the powers of the Euler operator and the constants, one gets, 
F{D) = {D + a){D + P), and the condition F{D)x''' = 0, gives, A = —a, —[3. The series 
solution is then 



lm=0 



-1 



d 



X- 



{D + a){D + l3) y^dx"^ ^ dx^ 



X 



(5) 



It is easily seen that, modulo an overall normalization factor, the above solution, is a rear 
ranged form of the well-known hypergeometric series. Making use of the following represen- 
tation of + /?), 

1 



Jo 



(6) 



and 



one can easily show that 

—1 d 



X 



1 



d"^ d , 
[x^^ + 7—) 



ml 



(8) 



{D + a) dx"^ ' ' dx 

Choosing the normahzation constant to match the conventional definition, we can write the 
normalized solution for the HGDE as, 



F(a.ft,;.) = (-l)-"^|^-|iJ2|.xp 



-1 I £ d'- 

I X h 7 

{D + a) \ dx^ dx ^ 



. X 



(9) 



The exponential form is a novel expression for F(q;, /3; 7; x), unknown in the literature, to 
the best of the authors' knowledge. The choice A = —a will lead again to the hypergcomctric 
series, as the series is symmetric under the exchange of a and (5. Since, j^^{xd^/dx'^ + 
jd/dx) lowers the degree of a given monomial by one, it can be seen that, when — /3 is an 
integer, the above is a polynomial solution of the HGDE. 

Similarly, the normalized polynomial solution for the confiuent hypergeometric (CH) 
differential equation. 



d^ , . d 



(10) 



IS 



$(q;;7;x) = ^ _ , exp 



rf2 ^ d' 
dx"^ ^ dx 



. X 



(11) 



r(7 - a) 

It is worth mentioning that, one can also modify the HGDE, by multiplying it with x, 
which yields, F{D) = D{D + j — l) and P{x,d/dx) = —x{D + a){D + P). The corresponding 
solution is 

-1 



F(a,/?;7;x) = qo,i-,) E(-l) 

Km=0 



-x{D + a){D + (3) 



X 



(0,1-7) 



, (12) 



_D(D + 7-l) 

since F{D)x'^ — gives A = 0, 1 — 7. For A = 0, the series can be written in the well-known 
form 



^0 (7)n n\ 



(13) 



where, 



a„ = a{a + l)(a + 2) • ■ ■ (a + — 1) 



(14) 



is the Pochammer symbol. The corresponding exponential form is 



F{a,f];'y;x) = Coexp 



X- 



. 1 



(15) 



(^ + 7) 

The other choice, namely A = 1— 7, yields the other linearly independent solution. Analogous 
results follow for CHDE. Hence, both the polynomial and the series solutions are obtained 
by the present method. 

The procedure for finding the solutions for CHDE and HGDE easily extends to the 
generalized hypergeometric cases0, given by the equations of the type. 



[6(6 + 61 - 1) • • • (6 + 6p - 1) - z(6 + ai) • • ■ (6 + Op + l)]y = , 



(16) 



where 6 = zd/dz, with z being, in general, complex. The generalized hypergeometric series, 
is denoted by p+i-Fp, where p + 1 and p are the number of parameters appearing in the 
numerator and the denominator of the generalized HG series, respectively. Since 6 is the 
Euler operator and the DE is already separated into F{Q) and P{z,d/dz), the solutions 
easily follow. We consider the case of the 3F2 series below, because of its importance in the 
quantum theory of angular momentum^. 
The DE for the 3F2 series 



[6(6 + 61 - 1)(6 + 62 - 1) - zie + ai)(6 + a2)(6 + ag)] 3F2 = 



(17) 



yields F(6) = 6(6 + 61 - 1)(6 + 62 - 1) and P{z, d/dz) 
The solution is then 



-z{Q + ai)(6 + a2)(6 + 03). 



3F2 - C(o,i-6i,i-;,2) < XI 



^ m=0 



1 



F(6) 



P{z,d/dz) 



z 



(0,1-61,1-62) 



(18) 



where F{Q)z'^ = 0, has yielded three solutions, A = 0, 1 — 61, 1 — 62- For A = 0, the above 
series can be expanded and rearranged in the conventional form 
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^ {ai)n{a2)n{a3)n Z"" , , 

^0 (&l)n(&2)n ri! 

Like the previous examples, an exponential form for the 3F2 can also be written down easily. 

The other roots, 1 — bi and 1 — &2, give the other two linearly independent solutions. It 
is clear from the expression for the 3F2 series, that it terminates, when either ai, 02 or 03, is 
a negative integer or zero. 

Multiplying the above DE with —1/z, one gets, 

[(0 + ai)(e + a2)(e + as) - ^(6 + 61 - l)(e + 62 - 1)] 3i^2 = , (20) 

where F{Q) and P{z,d/dz) are now given by F(6) = (6 + ai)(0 + a2)(0 + 03) and 
P{z,d/dz) = —{d/dz){Q + 61 — 1)(0 + &2 — !)• Therefore, the solution, in this case, turns 
out to be 

1 



3F2 - C(ai,a2,a3) { X! (^^^ 
lm=0 



-P{z,d/dz) 



where ai, 02, 03 have to be negative integers to yield a polynomial solution. 

By now, it is clear that for a number of DEs, F{D)x^ = 0, leads to the linearly indepen- 
dent solutions, when the solutions are nondegenerate. In cases, where this is not possible, 
appropriate modification of the DE, has straightforwardly yielded the linearly independent 
solutions. Below, we give two more examples to illustrate these points. In particular, the 
second example deals with the case, where F{D)x^ = 0, leads to degenerate solutions. These 
examples will also point out the connection of the indicial equation, in the conventional se- 
ries solution method, to the present one. In the standard series solution methodEi, for the 
DE 

, d^V dy , , 

Ax-^ + 2-^ + y = , (22) 

the two distinct roots of the indicial equation c = 0, 1/2, lead to two linearly independent 
solutions. Multiplication by x yields 

Ax'p^ + 2x^ + xy = , (23) 
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which imphes, F{D) = Ax'^-^+2xj^ = 2D{2D-1) and P{x, d/dx) = x. Hence, F{D)x^ = 
obtains A = 0, 1/2, which provide the two hnearly independent solutions, as is obtained in 
the standard approach. 
In our second example, 

the roots of the indicial equation, = are degenerate. Multiplying the above DE with x, 
to bring it to the form given by Eq. one gets (D^ + x)y = 0; hence F{D) = and 
P{x) = X. It is clear that, F{D)x^ = 0, also leads to the same degenerate case as obtained 
by the method of series solution. In this scenario, one has to employ the established methods 
for finding out the other linearly independent solutionlll. 

The procedure developed here is applicable to a wide range of functions and polynomials. 
Some of the well-known ones, explicitly checked by the authors are, Meijer's G-function, 
Struve, Lomel, Anger, Weber, Bessel functions, Gegenbauer, Neumann's, Jacobi, Schlafii, 
Whittaker and Chebyshev polynomials. It is also applicable in the periodic cases. For 
example, solution of DE with the following periodic potential 

d^y , . 

— ■ + a cos{x)y = , (25) 
dx'^ 

can be found, after multiplying Eq. (^) by and rewriting x"^-^ as {D — 1)D, 



\ {2n,)\ J \^^^{2[m + X/2 + l-r + Y.T=i 

Xnrn+j:Z,n,+X/2) (26) 



Here, A = or 1. In the same manner, one can write down the solutions for the Mathieu's 
equation as welll. 

Although this paper is devoted to the single variable case, it is worth pointing out that 
the generalization of this method to the many-variable case is immediate^'i. This can be 
accomplished by denoting D = Y^iDi = Y.iXi-£r, where i = 1,2,- --N. Using the fact. 
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F{D)X'^ = has solutions, in the space of monomial symmetric functional^!, a number of 
many-body equations can be solved, in a manner analogous to the single variable case. 



We now proceed to study the algebraic properties of the space of solutions. The ad- 
vantages of the present approach, as compared to the previous oneJi, lie in the following 
two facts. First of all, a priori, no symmetry of the DE is assumed. Secondly, the ladder 
operators are straightforward to construct in the space of the monomials, which can be 
brought to the space of solutions, via a similarity transformation, with the aid of the expo- 
nential form of the solutions. The only criterion in the choice of the ladder operators, in the 
space of monomials, is that, after the similarity transformation, the resulting operators, are 
well-defined; these operators then yield the symmetry algebra. 

Keeping in mind, the appearance of CHDE and HGDE in diverse physical systems, 
we explicitly work out the generators of the symmetry algebras, for these two cases, for 
the corresponding polynomial solutions. Application to other DEs can be carried out in a 
similar manner. 



leads to a polynomial solution, only if, a is a negative integer {—n) or zero. The above form 
of the solution, immediately suggests the lowering operator to be. 



since J_ commutes with the exponential and hence can lead to a lowering operator at the 
level of the polynomial. J_ reduces the degree of the polynomial by one. 



III. ALGEBRAIC STRUCTURE OF THE SOLUTION SPACE 



The solution of the CHDE, 




(27) 




(28) 



J_e = n(7 + n - l)e x 



,n-l 



(29) 



After taking into account the normalization factors, we get 
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Choosing J+ = x, as the raising operator at the level of the monomials and introducing an 
identity operator, in the following manner, 

e-^-xx" = e-^-xe+'^-e-^-x'' , (31) 

one gets, 

$(-n;7;a;) = -(n + 7)<l>(-n - l;7;a;) . (32) 



^ d rf2 ^■ 

X - 2x- 7 + x—^ + 7— 

ax ax^ ax 



The similarity transformation of x, the raising operator, at the level of the monomials, 
led to the raising operator at the level of the polynomials. It can be easily seen that, 
[J+, J_] = —2 Jo, where Jo — D + 7/2 and 

[Jo,J±]=±J± . (33) 

This is the well-known SU{1,1) dynamical algebra, in the solution space of the CHDE. 
With an appropriate choice of the measure, the operators J+ and J_ can be made the 
formal adjoints of each other. It is worth pointing out that, the simplest choice d/dx as the 
lowering operator for the monomials, does not lead to a convergent expression at the level 
of the polynomials. It should be pointed out that other forms of ladder operators can also 
be found. 

A straightforward calculation, taking x and d/dx as the raising and lowering operators 
at the level of the monomials, leads to the Heisenberg algebra, [a, a^] = 1, for the Hermite 
polynomials 

//„(x) = C„exp(-i^).x" . (34) 

The exponential form of the solution of HGDE suggests the simplest lowering operator to 
be 

1 / d^ d\ ~ 
10 



which lowers the degree of the monomials by one. Here T = l/{D+a) and J_ = [x{(P / dx'^) + 
'-j{d/ dx)]. It is to be noted that J_ can also act as the lowering operator at the level of the 
monomials; however, the latter choice will not lead to a convergent expression in the solution 
space. 

Operating J_ on x", one gets, 

J_x"= "/^ + "-V x"-^ , (36) 
(a + n-1) ^ ^ 

which straightforwardly extends to, 

1 ( d^ d\ 

X— ^ + 7— F(q;, -n; 7; x) = -nF(a;; -n + 1; 7; x) , (37) 



{D + a) \ dx"^ dx ^ 

at the level of the polynomials. The raising operator, at the level of the monomials, needs 
to be chosen carefully due to the presence of J_ in the exponential. 

At this point, we note that, the existence of a canonical conjugate operator for J__, 
such that, [J_, J+] = 1, akin to the Heisenberg algebra, would imply 

e-^'- J+x" = (J+ - l)e-^"-x" . (38) 

This suggests J+ can be used as a raising operator in the space of monomials, which can be 
written in a compact form in the space of solutions of the HGDE. 

The procedure for finding the canonical conjugate is straightforward^. We briefly outline 
the same below and use it for finding J+. Although we have used this procedure for conve- 
nience, other forms of raising operators are possible to construct. The method for obtaining 
them will be illustrated below. Denoting Jq = xdj dx^'y jl^ one finds [Jo, J-] = — J-. Start- 
ing with a function T{Jq), whose required properties will become clear in the subsequent 
steps, we define 

J+ = xT(Jo) 

= J+T{Jo) . (39) 



It is easy to check that. 
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[Jo, J+] = J+ and [J+, J- 



-2Jn 



(40) 



Since — 2Jo = gi^Jo) —g{Jo — l), with g{Jo) = —Jo{Jo + l), the Casimir operator C commuting 
with all the three generators can be written as, C = J_J^ + g{Jo) = J+J- + g{Jo — 1)0. 
Starting from [J_, J^] = 1, one finds 

f{Jo)T{Jo)J-J+ - T{Jo - l)f (Jo - 1) J+ J- = 1 , 



which leads to 



T(Jo) = [f (Jo)] 



_i Jo + S 



[C-g{Jo)] ■ 
Here 6 is an arbitrary constant. Simplification yields 

~ ^ ( Jo + a - 1) (Jo + ^ - 1) 

^ (Jo + 7-l)^o ■ 

Demanding that [J_, J+] = 1, holds on the lowest monomial i.e., 1, one obtains 5 = 1. This 
leads to 



^ (Jo + g-l) ^^ 



(41) 



(^0 + 7-1) 

Now the raising operator in the solution space can be obtained via a similarity transforma- 
tion: 



e-'^-J+e'^-e-'^-x" 



[a + n) 



(7 + n) 

The above expression after restoring the normalization factors yields. 



(42) 



\ _ (^0 + ^-1) ^' 
(^0 + 7-1) . 



F(a; — n; 7; x) = F{a; —n — 1; 7; x) 



(43) 



As noted earlier, the raising operator obtained above is not unique, one can construct other 
raising operators; some examples are 



J+ = (x + x^—) and J+ = [x + x^ — )T{C, Jq) 
ax ax 



For the former case, explicit computation leads to 



(44) 
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1 - -—{x + X — ) 

D + dx' 



F{a; —n; 7; x) = F{a; —n — 1; 7; x) 



(45) 



It is interesting to note that, the algebra satisfied by J+, J_ and D is a. quadratic algebr 



since 



[J+, J_ 



-2(7 + 1/2)D - - 7 



(46) 



and 



±J. 



±- 



The algebra satisfied by J+, J_ and Jo is the well known SU{1, 1) algebra. Here J+, J_ can 
be made formal adjoints of each other with the appropriate choice of the measure. The point 
to note is that, both the above mentioned SU{1, 1) and the quadratically deformed algebra, 
have the solution space of HGDE, as their irreducible representations. 

The above ladder operators can be suitably rearranged to yield the ladder operators 
obtained by the factorization method (FM)0'ii and SUSY-QM0. These ladder operators 
be used for obtaining wavefunctions of a number of quantum mechanical problems&B, 



can 



after appropriate measures are introduced. We take the well-known example of the Laguerre 
DE, for the purpose of establishing the above mentioned connection: 



x-r^ + (a + 1 - xj— + n 
dx^ dx 







(47) 



The solution of the above DE is the familiar Laguerre polynomial, 

(-ir 



exp 



d^ ( + 1) 
dx'^ dx 



X 



(48) 



In order to find the raising operator for L^, analogous to the ones obtained by FM and 
SUSY-QM, we start with 



d 

X + X- n 

dx 



as the raising operator at the level of the monomials. Following the method employed for 
CHDE, we get 
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d . , 
X- — \- [a + 1) + n — X 
ax 



L: = {n + l)L 



n+l 



(49) 



This can be cast in the form 



n 

k=l 



X- — \- a + k — X 
ax 



(50) 



or 



LI = A\a + l)A\a + 2) ■ ■ ■ A\a + n)l 



where A^{a + 



n] 



x-^ + a + n — X 



. These shifted operators have found apphcation in 
the construction of coherent states^. This form of the polynomial, modulo normalizations; 
matches with the ones obtained from FM and SUSY-QM. Using the raising operator for 
CHDE and HGDE, one can al;so obtain similar expressions for other polynomials. 

It should be emphasized that FM and SUSY-QM, have taken recourse to a special prop- 
erty of the equations under study called shape invariance@ for arriving at these results. 
Furthermore it was necessary to introduce a host of intermediate systems. In contrast the 
present technique does not presume any such special properties of the DE. 

Since we have already provided the raising operator, for the sake of completeness, we 
derive the lowering operator for in a convenient form. Starting from. 



d d"^ , d 

X- n + x—^ + (a + 1) — 

dx dx"^ dx 



a suitable similarity transformation yields 



x"' = n(n + a)x"- ^ 



d 



X- 



dx 



n 



l = -{n + a)LU ■ (51) 

Eq. ( ^91) and Eq. (|5T|) are the standard recurrence relations! for the Laguerre polynomials. 
We can also derive those operators which change the values of a. The following steps lead 
to, 

-1)" d 



n\ dx 



exp 



d"^ ( + 1) 
dx"^ dx 



X 



(52) 
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i-iy 



ni 



exp 



dx dx 



d 

dx 



X" 



--KTlix) . 



(53) 



Similarly for HG series we get, 



-^F{a,(3;-f;x) 



X exp 



, r(a-/3)r(7) 
'^^ ^ r(7-/3)r(a) 
-1 



X 



_( d^ , dT" 
{D + a + 1) ydx^ ^^^^ 'dx, 



. X 



(54) 



= — F(q; + 1,/3 + 1;7 + 1;x) 
7 

IV. APPPLICATION TO CLASSICAL ORTHOGONAL POLYNOMIALS 



(55) 



A. Rodriguez Formula 



In this section, wc elaborate on tiie applicability of the approach developed here for 
finding other properties of the solution space. We start with the Rodriguez formula (RF) and 
show how these can be obtained, with the examples of Laguerre and Hermite polynomials. 
RF for Laguerre polynomials 



For simplicity, we consider the expression of L°, 

rO/ ^ (-1)" ( d\ , 

Defining B = xd^/dx"^ + d/dx and introducing an identity operator, we get 



(56) 



-1)" 



The above expression simplifies to 

f-1 



Llix) 



d d^ d 

X - 2x- 1 + x-p-r + — 

dx dx'^ dx 



(57) 



which can be written in the form 
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Further simplification yields, 



d d 

dx dx 



and hence 



(58) 



1 r/" 



(59) 



This is the well-known Rodriguez formula for the Laguerre polynomials. 
RF for the Hermite polynomials 

The normalized solution of the Hermite DE is given by, 



i/„(x) = 2"exp(-i^)x" . 



Denoting A = and writing 



(60) 



HJx) = 2"e-^xVe-^l 



(61) 



we obtain 



Hnix) 



d ' 

2x- — 

dx 



Introducing an identity operator in the form. 



HJx) = e^\-^" 



d 

2a; - — 

dx 



one obtains the RF as 



(62) 



(63) 



2 d?' _„2 



Hn{x) = -l^e^ — e 
dx'^ 



(64) 



Eq. (|57D and Eq. ( |62D reveal that simplification of the raising operators by introducing 
appropriate identity operators led to the RF. Using the raising operators obtained in the 
previous section, one can easily extend these results to other polynomials. 



16 



One can also start from the RF and obtain Eg. (|60|) . Explicitly, 



Hjx) 



-IT 



d 

— -2x 

ax 



_ 1 d 
^ ~ 2dx 



e-^e^ 1 



(65) 



The above expression, by suitable manipulations, leads to 

H^{x) = 2-eM-\^)^'' ■ 

This shows the procedure to obtain the exponential form of the solutions, starting from 
the known RF for a given solution. 



B. Generating Functions 

Below, we outline the method of getting the generating functions (GF) taking Laguerre 
and Chebyshev as examples. The advantage of having an exponential form for the solution 
comes out naturally through these examples. 
GF for the Laguerre Polynomial. 

Defining the generating function as 



g{x,t) = Y.Ll{x)e 



n=0 



and substituting the expression for L^{x) in the above equation, we get, 

oo { -\\n 



n=0 



n! 



Action i? on e ^* is easy to compute: 
g{x,t) = e"^'* 



-xt 



1 - {xt^ -t) + ^{xH^ - 4xt^ + 2t^) 



(66) 



(67) 



(68) 
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This series can be summed and written in the compact form 

_ exp[-xt/(l-t)] 

g{x,t) is the well-known GF for L'^ix). The procedure outlined above can be straightfor- 
wardly extended to include, a ^ cases also. 
Chebyshev type II 

We consider type II Chebyshev polynomial for illustration. The solution of Chebyshev 
type II DE can be written as 

1 £ 



Unix) = 2"exp 
and its GF is defined as, 



2{D + n + 2) dx^ 



X 



g{x,t) = J2Unix)f^ 

n=0 



Substituting the expression for Un{x) in the above equation, we obtain, 

1 ^2 

n=0 



gix,t) = J22^exp 



' 2{D + n + 2) dx"^ 



X t 



(69) 



(70) 



l + 2xt + (4x^ - l)t2 + 



(71) 



which can be summed to yield the GF as 

g{x,t) = 



1 - 2xt + t2 

This procedure for obtaining GF can be generalized to other orthogonal polynomials. It is 
clear that the exponential form of the solutions enables one to find the GFs straightforwardly. 

Generalized GFs find application in the construction of coherent and squeezed stateslllS. 
It is worth mentioning that the exponential form of the Hermite polynomials has been 
connected with the Gauss transform^. Our results provide analytic expressions for the 
Gauss transform, with appropriate parameter value, of various polynomials. Construction 
of generalized coherent states, associated with these polynomials is currently under progress 
and will be reported elsewhere. 



V. CONCLUSIONS. 



In conclusion, the solution of a wide class of liner differential equations, which can be 
cast into a form, where a function of the Euler operator and constants separates from the 
rest, can be written in a closed form, which makes the algebraic properties of the solution 
space transparent. Explicit examples dealing with confluent hypergeometric, hypergeometric 
and generalized hypergeometric equations were analyzed, where, not only novel expressions 
connecting the solution space, with the space of the monomials were written down, but also, 
utilized for unravelling the dynamical symmetryies underlying the solution spaces. 

The fact that, a priori no assumption was made about the symmetry of the equation 
understudy, makes this approach attractive. Although, we have analyzed here well-known 
examples, an exhaustive study reveals that the present approach extends to a host of other 
functions and polynomials. Some of these are, Meijer's G-Function, Struve, Lomel, Anger, 
Weber, Bessel functions, Gegenbauer, Neumann's, Jacobi, Schlafli, Whittaker, Chebyshev 
polynomials. These functions and polynomials manifest in diverse branches of physics and 
mathematics. The novel expressions for the solutions presented here will help in unravelling 
various properties of these functions and polynomials. 

Here, we have only briefly mentioned about the multivariate cases dealing with correlated 
systems. The solution space of these equations have rich symmetry and they are connected 
with random matrices which find application in diverse areasil. Hence, a deeper analysis of 
the the problem is also warranted. Furthermore, the connection of this approach with the 
Gauss transform can be used for constructing generalized coherent states. It is worth men- 
tioning that, the exponential form of the Hermite polynomial has already found application 
in the construction of coherent and squeezed statesSS. 

Apart from finding exact solutions, the present approach can also be utilized for finding 
approximate solutions to differential equationsil. Some of these works are currently under 
progress and we hope to report the findings in the near future. 
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